The role of mutation rate variation and genetic diversity in the architecture of human disease

Background We have investigated the role that the mutation rate and the structure of genetic variation at a locus play in determining whether a gene is involved in disease. We predict that the mutation rate and its genetic diversity should be higher in genes associated with disease, unless all genes that could cause disease have already been identified. Results Consistent with our predictions we find that genes associated with Mendelian and complex disease are substantially longer than non-disease genes. However, we find that both Mendelian and complex disease genes are found in regions of the genome with relatively low mutation rates, as inferred from intron divergence between humans and chimpanzees, and they are predicted to have similar rates of non-synonymous mutation as other genes. Finally, we find that disease genes are in regions of significantly elevated genetic diversity, even when variation in the rate of mutation is controlled for. The effect is small nevertheless. Conclusions Our results suggest that gene length contributes to whether a gene is associated with disease. However, the mutation rate and the genetic architecture of the locus appear to play only a minor role in determining whether a gene is associated with disease

Evidence for Pervasive Adaptive Protein Evolution in Wild Mice

The relative contributions of neutral and adaptive substitutions to molecular evolution has been one of the most controversial issues in evolutionary biology for more than 40 years. The analysis of within-species nucleotide polymorphism and between-species divergence data supports a widespread role for adaptive protein evolution in certain taxa. For example, estimates of the proportion of adaptive amino acid substitutions (alpha) are 50% or more in enteric bacteria and Drosophila. In contrast, recent estimates of alpha for hominids have been at most 13%. Here, we estimate alpha for protein sequences of murid rodents based on nucleotide polymorphism data from multiple genes in a population of the house mouse subspecies Mus musculus castaneus, which inhabits the ancestral range of the Mus species complex and nucleotide divergence between M. m. castaneus and M. famulus or the rat. We estimate that 57% of amino acid substitutions in murids have been driven by positive selection. Hominids, therefore, are exceptional in having low apparent levels of adaptive protein evolution. The high frequency of adaptive amino acid substitutions in wild mice is consistent with their large effective population size, leading to effective natural selection at the molecular level. Effective natural selection also manifests itself as a paucity of effectively neutral nonsynonymous mutations in M. m. castaneus compared to humans

Identifying the favored mutation in a positive selective sweep.

Most approaches that capture signatures of selective sweeps in population genomics data do not identify the specific mutation favored by selection. We present iSAFE (for "integrated selection of allele favored by evolution"), a method that enables researchers to accurately pinpoint the favored mutation in a large region (∼5 Mbp) by using a statistic derived solely from population genetics signals. iSAFE does not require knowledge of demography, the phenotype under selection, or functional annotations of mutations

Molecular Characterization of the Gastrula in the Turtle Emys orbicularis: An Evolutionary Perspective on Gastrulation

Due to the presence of a blastopore as in amphibians, the turtle has been suggested to exemplify a transition form from an amphibian- to an avian-type gastrulation pattern. In order to test this hypothesis and gain insight into the emergence of the unique characteristics of amniotes during gastrulation, we have performed the first molecular characterization of the gastrula in a reptile, the turtle Emys orbicularis. The study of Brachyury, Lim1, Otx2 and Otx5 expression patterns points to a highly conserved dynamic of expression with amniote model organisms and makes it possible to identify the site of mesoderm internalization, which is a long-standing issue in reptiles. Analysis of Brachyury expression also highlights the presence of two distinct phases, less easily recognizable in model organisms and respectively characterized by an early ring-shaped and a later bilateral symmetrical territory. Systematic comparisons with tetrapod model organisms lead to new insights into the relationships of the blastopore/blastoporal plate system shared by all reptiles, with the blastopore of amphibians and the primitive streak of birds and mammals. The biphasic Brachyury expression pattern is also consistent with recent models of emergence of bilateral symmetry, which raises the question of its evolutionary significance

Widespread Genomic Signatures of Natural Selection in Hominid Evolution

Selection acting on genomic functional elements can be detected by its indirect effects on population diversity at linked neutral sites. To illuminate the selective forces that shaped hominid evolution, we analyzed the genomic distributions of human polymorphisms and sequence differences among five primate species relative to the locations of conserved sequence features. Neutral sequence diversity in human and ancestral hominid populations is substantially reduced near such features, resulting in a surprisingly large genome average diversity reduction due to selection of 19–26% on the autosomes and 12–40% on the X chromosome. The overall trends are broadly consistent with “background selection” or hitchhiking in ancestral populations acting to remove deleterious variants. Average selection is much stronger on exonic (both protein-coding and untranslated) conserved features than non-exonic features. Long term selection, rather than complex speciation scenarios, explains the large intragenomic variation in human/chimpanzee divergence. Our analyses reveal a dominant role for selection in shaping genomic diversity and divergence patterns, clarify hominid evolution, and provide a baseline for investigating specific selective events

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

MM received funding from the Max Planck innovation funds for this project. PGDF was supported by a Marie Curie European Reintegration Grant (proposal nr 270891). CE was supported by German Science Foundation grants (DFG, EI 841/4-1 and EI 841/6-1)

A multilocus assay reveals high nucleotide diversity and limited differentiation among Scandinavian willow grouse (Lagopus lagopus)

<p>Abstract</p> <p>Background</p> <p>There is so far very little data on autosomal nucleotide diversity in birds, except for data from the domesticated chicken and some passerines species. Estimates of nucleotide diversity reported so far in birds have been high (~10^-3) and a likely explanation for this is the generally higher effective population sizes compared to mammals. In this study, the level of nucleotide diversity has been examined in the willow grouse, a non-domesticated bird species from the order Galliformes, which also holds the chicken. The willow grouse (<it>Lagopus lagopus</it>) has an almost circumpolar distribution but is absent from Greenland and the north Atlantic islands. It primarily inhabits tundra, forest edge habitats and sub-alpine vegetation. Willow grouse are hunted throughout its range, and regionally it is a game bird of great cultural and economical importance.</p> <p>Results</p> <p>We sequenced 18 autosomal protein coding loci from approximately 15–18 individuals per population. We found a total of 127 SNP's, which corresponds to 1 SNP every 51 bp. 26 SNP's were amino acid replacement substitutions. Total nucleotide diversity (<it>π</it>_<it>t</it>) was between 1.30 × 10^-4and 7.66 × 10^-3(average <it>π</it>_{<it>t </it>}= 2.72 × 10^-3± 2.06 × 10^-3) and silent nucleotide diversity varied between 4.20 × 10^-4and 2.76 × 10^-2(average <it>π</it>_{<it>S </it>}= 9.22 × 10^-3± 7.43 × 10^-4). The synonymous diversity is approximately 20 times higher than in humans and two times higher than in chicken. Non-synonymous diversity was on average 18 times lower than the synonymous diversity and varied between 0 and 4.90 × 10^-3(average <it>π</it>_{<it>a </it>}= 5.08 × 10^-4± 7.43 × 10³), which suggest that purifying selection is strong in these genes. <it>F</it>_STvalues based on synonymous SNP's varied between -5.60 × 10^-4and 0.20 among loci and revealed low levels of differentiation among the four localities, with an overall value of <it>F</it>_ST= 0.03 (95% CI: 0.006 – 0.057) over 60 unlinked loci. Non-synonymous SNP's gave similar results. Low levels of linkage disequilibrium were observed within genes, with an average r²= 0.084 ± 0.110, which is expected for a large outbred population with no population differentiation. The mean per site per generation recombination parameter (ρ) was comparably high (0.028 ± 0.018), indicating high recombination rates in these genes.</p> <p>Conclusion</p> <p>We found unusually high levels of nucleotide diversity in the Scandinavian willow grouse as well as very little population structure among localities with up to 1647 km distance. There are also low levels of linkage disequilibrium within the genes and the population recombination rate is high, which is indicative of an old panmictic population, where recombination has had time to break up any haplotype blocks. The non-synonymous nucleotide diversity is low compared with the silent, which is in agreement with effective purifying selection, possibly due to the large effective population size.</p

Recurrent Modification of a Conserved Cis-Regulatory Element Underlies Fruit Fly Pigmentation Diversity

The development of morphological traits occurs through the collective action of networks of genes connected at the level of gene expression. As any node in a network may be a target of evolutionary change, the recurrent targeting of the same node would indicate that the path of evolution is biased for the relevant trait and network. Although examples of parallel evolution have implicated recurrent modification of the same gene and cis-regulatory element (CRE), little is known about the mutational and molecular paths of parallel CRE evolution. In Drosophila melanogaster fruit flies, the Bric-à-brac (Bab) transcription factors control the development of a suite of sexually dimorphic traits on the posterior abdomen. Female-specific Bab expression is regulated by the dimorphic element, a CRE that possesses direct inputs from body plan (ABD-B) and sex-determination (DSX) transcription factors. Here, we find that the recurrent evolutionary modification of this CRE underlies both intraspecific and interspecific variation in female pigmentation in the melanogaster species group. By reconstructing the sequence and regulatory activity of the ancestral Drosophila melanogaster dimorphic element, we demonstrate that a handful of mutations were sufficient to create independent CRE alleles with differing activities. Moreover, intraspecific and interspecific dimorphic element evolution proceeded with little to no alterations to the known body plan and sex-determination regulatory linkages. Collectively, our findings represent an example where the paths of evolution appear biased to a specific CRE, and drastic changes in function were accompanied by deep conservation of key regulatory linkages. © 2013 Rogers et al

Genomic and expression analysis of multiple Sry loci from a single Rattus norvegicus Y chromosome

BACKGROUND: Sry is a gene known to be essential for testis determination but is also transcribed in adult male tissues. The laboratory rat, Rattus norvegicus, has multiple Y chromosome copies of Sry while most mammals have only a single copy. DNA sequence comparisons with other rodents with multiple Sry copies are inconsistent in divergence patterns and functionality of the multiple copies. To address hypotheses of divergence, gene conversion and functional constraints, we sequenced Sry loci from a single R. norvegicus Y chromosome from the Spontaneously Hypertensive Rat strain (SHR) and analyzed DNA sequences for homology among copies. Next, to determine whether all copies of Sry are expressed, we developed a modification of the fluorescent marked capillary electrophoresis method to generate three different sized amplification products to identify Sry copies. We applied this fragment analysis method to both genomic DNA and cDNA prepared from mRNA from testis and adrenal gland of adult male rats. RESULTS: Y chromosome fragments were amplified and sequenced using primers that included the entire Sry coding region and flanking sequences. The analysis of these sequences identified six Sry loci on the Y chromosome. These are paralogous copies consistent with a single phylogeny and the divergence between any two copies is less than 2%. All copies have a conserved reading frame and amino acid sequence consistent with function. Fragment analysis of genomic DNA showed close approximations of experimental with predicted values, validating the use of this method to identify proportions of each copy. Using the fragment analysis procedure with cDNA samples showed the Sry copies expressed were significantly different from the genomic distribution (testis p < 0.001, adrenal gland p < 0.001), and the testis and adrenal copy distribution in the transcripts were also significantly different from each other (p < 0.001). Total Sry transcript expression, analyzed by real-time PCR, showed significantly higher levels of Sry in testis than adrenal gland (p, 0.001). CONCLUSION: The SHR Y chromosome contains at least 6 full length copies of the Sry gene. These copies have a conserved coding region and conserved amino acid sequence. The pattern of divergence is not consistent with gene conversion as the mechanism for this conservation. Expression studies show multiple copies expressed in the adult male testis and adrenal glands, with tissue specific differences in expression patterns. Both the DNA sequence analysis and RNA transcript expression analysis are consistent with more than one copy having function and selection preventing divergence although we have no functional evidence